A cosmological concordance model with dynamical vacuum term 
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We demonstrate that creation of dark-matter particles at a constant rate implies the existence of 
a cosmological term that decays linearly with the Hubble rate. We discuss the cosmological model 
that arises in this context and test it against observations of the first acoustic peak in the cosmic 
microwave background (CMB) anisotropy spectrum, the Hubble diagram for supernovas of type la 
(SNIa), the distance scale of baryonic acoustic oscillations (BAO) and the distribution of large scale 
structures (LSS). We show that a good concordance is obtained, albeit with a higher value of the 
present matter abundance than in the ACDM model. We also comment on general features of the 
CMB anisotropy spectrum and on the cosmic coincidence problem. 

I. INTRODUCTION 

The cosmological constant problem is usually described as a huge difference between the vacuum energy density 
derived by quantum field theories and the observed value of the cosmological constant. The energy density associated 
to the vacuum fluctuations of free fields is a divergent quantity, and any natural cutoff we use to regularize it (for 
example, the Planck mass or the energy scale of the quantum chromodynamic (QCD) phase transition, the latest 
cosmological vacuum transition) leads to values many orders of magnitude above the observed A. But the applicability 
of flat-spacetimc quantum field theory in the expanding universe is doubtful, since a non- vanishing bare cosmological 
constant is not compatible with a flat spacetime. The vacuum contribution to (semi-classical) gravity is only properly 
calculated when one considers a non-flat background and renormalizes the divergent contributions in some appropriate 
way. This is a difficult task and, in general, may depend on the renormalization method used. In the simpler case of 
free conformal fields in a de Sitter spacetime [l[ , the renormalized vacuum density is proportional to H 4 , where H is 
the Hubble parameter. In the high-energy limit this result can be used to obtain a non-singular cosmology, with an 
initial quasi-de Sitter phase, giving origin to a subsequent radiation dominate period @, Q. However, in the present 
universe the above result leads to a very tiny vacuum term. On the other hand, it may well be possible that the 
observed cosmological term is just a cosmological constant of purely geometric nature. In this case there would be no 
natural (micro-) physical scale associated to it. Furthermore, both the gravitational action of the quantum vacuum 
and the approximate coincidence between such a constant and the present matter density (the cosmic coincidence 
problem) would remain unexplained. 

As long as this issue cannot be solved on a fundamental level, (semi-) phenomenological models which avoid the 
theoretical difficulties pointed out above may help to get further insight into the problem. In this letter we argue 
that, by modeling the cosmological constant in terms of a decaying vacuum energy, it is possible to obtain a model 
which is competitive with the standard one from the observational point of view. A crucial point here is to consider 
the contribution of interacting fields to the vacuum energy. Indeed, some authors have been arguing that the energy 
density associated to the QCD vacuum condensate in a low-energy de Sitter (or approximately de Sitter) spacetime 
is given by 

A w m 3 H, (1) 

where m rs 150 MeV is the energy scale of the QCD phase transition Q. In the present universe this density 
is dominant when compared to the contribution of free fields, if 4 , and has the correct order of magnitude of the 
observed cosmological term. In this paper we will trace back this scaling law A oc H to a process of matter production 
in the expanding spacetime. 

II. PARTICLE PRODUCTION AND THE VACUUM DENSITY 

Let us consider a low-energy, spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime, where non- 
relativistic dark particles of mass M are created out of the vacuum at a constant rate T. The particle number balance 
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(sometimes called Boltzman equation) for this process will be given by 

^ t (a*n)=rn, (2) 
where n is the particle number density at a given time. It can be rewritten as 

p m + 2>Hp m = Tp m , (3) 
where p m — nM is the matter density. This equation is characteristic of models with interaction in the dark sector 

This particle production from vacuum fluctuations necessarily involves a gravitational backreaction. In addition to 
equation ([3]) we have the Friedmann equation 

p = Pm +A = 3H 2 , (4) 

where A is the vacuum contribution to the energy density. We shall assume p\ = —A for the vacuum pressure, 
admitting, however, A to be time varying. The total energy density satisfies the conservation equation 

p + 3H(p + p)=0, (5) 

provided we take A = 2TH + A , where Ao is a constant of integration, which corresponds to a purely geometric cos- 
mological constant. Since there is no natural (microphysical) energy scale associated to this constant, it is disregarded 
here. In this way we have 

A = 2YH, (6) 

and equation ([5]) can be rewritten as 

p m + 3Hp m = -A. (7) 

We can see that the particle creation is concomitant with a decay in the vacuum energy density. Equation ([6]) is 
the same time variation law ([IJ obtained in a quasi-de Sitter spacctime when the contribution of the QCD condensate 
is taken into account [![. In this case we have T sa to 3 , where m is the energy scale of the QCD phase transition 1 . As 
discussed in Q , the association of the QCD energy scale with the observed vacuum energy is natural and, in some way, 
expected, since the QCD confinement, with breaking of chiral symmetry, is the last cosmological vacuum transition we 
know. The QCD condensate formed after this low energy transition gives, in a de Sitter (or approximately de Sitter) 
spacctime, the correct value of the cosmological term (see below). Furthermore, owing to the maximal symmetry of 
the de Sitter background, the condensate pressure is the negative of its energy density, i.e., the same equation of state 
associated to a cosmological constant. 

Dividing equation ([6]) by 3H 2 , using K/(3H 2 ) = 1 — fl m and specifying to the present time (subscripts 0), we obtain 

r = \ (1 - Slmo) H . (8) 

In a de Sitter universe (fi m = 0) we would have T = 3H/2. That is, the particle creation rate would be equal (apart 
from a numerical factor) to the Gibbons-Hawking temperature associated to the event horizon |{|. On the other hand, 
by taking f2 mu « 1/3 for our present universe, we have H « m 3 and A « m 6 . The former relation is an expression 
of the old Eddington-Dirac large number coincidence [l(| , while the latter (known as Zeldovich's relation [ll| ) gives 
the correct order of magnitude for A. The exact value of T in © (and hence of Ao) will be given by observations in 
the next sections. Once T is fixed, no further fine-tuning will be necessary. The cosmological model to be built will 
have the same number of free parameters as the ACDM model, namely £l m o and Hq. 



III. THE MODEL 



From this section on we will analyze the cosmological model arising from a vacuum energy density scaling with 
time in accordance with the ansatz discussed above. When compared with the most precise observational constraints, 
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the model will show a very good concordance and, in some cases, less tensions than in the standard case. Therefore, 
besides providing a possible physical basis for the cosmological term, in the lines of the previous sections, it will also 
show to be competitive from the observational point of view. 

With A = 2TH we obtain, from the Friedmann equations, the solution [12, HH 



1/2 

i - n m0 + Ono(i + zf /2 +n r0 (i + z) 4 \ , (9) 



H 

Ho " 

where £l m o is the present relative matter density, and we have added conserved radiation with present density param- 
eter Q r Q. As discussed in (l3l - [l5| . for non-zero f2 r o the expression ^ is an approximate solution, differing only 1% 
from the exact one, since il r o ~ 8 x 10~ 5 <C 1. For Q r o = 0, the solution (j9)) is exact. 

For early times we obtain H 2 (z) = H 2 n r oz 4 , and the radiation era is indistinguishable from the standard one. 
On the other hand, for high redshifts the matter density scales as p m {z) = SHqQ^qZ 3 . The extra factor 17 m0 - as 
compared to the ACDM model - is owing to the late-time process of matter production. In order to have nowadays 
the same amount of matter, we need less matter in the past. Or, in other words, if we have the same amount of 
matter in the past (say, at the time of matter-radiation equality), this will lead to more matter today. We can also 
see from © that, in the asymptotic limit z —± — 1, the solution tends to the de Sitter solution. 

Note that, like the ACDM model, the above model has only two free parameters, namely £l m o an d i?o- On the 
other hand, it can not be reduced to the ACDM case except for z — > — 1. In this sense, it is falsifiable, that is, it may 
be ruled out by observations. 



IV. BACKGROUND TESTS 



The Hubble function (JSJ is all we need to test the model against observations of SNIa and BAO [13l4l5| . However, 
when testing the position of the first peak of the CMB anisotropy spectrum, some care is needed. The relation between 
the observed position of the peak and the acoustic scale I a is given by [l6| 



r 



h = l A (l-5 1 ), where 5 l = 0.267 [—) , (10) 



o.i 



with r = Pr/pm evaluated at the redshift of last scattering zi s . The acoustic scale is defined, as usual, by 

dz , f°° c s dz 







with the sound velocity 
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-1/2 

(12) 



where flbo and J7 7 o stand for the present density parameters of baryons and photons, respectively. The relation (|10[) 
does not depend on the dark energy model, but only on pre-recombination physics. However, the ratio r is, in our 
case, given by 

r = Tyj-Zls, (13) 
S2 m0 

with the extra factor fi m0 as compared to the ACDM expression. Note that, in all the above expressions, there is a 
dependence on Hq through the relative densities of baryons, photons and radiation, as well as through the redshift of 
last scattering. In our analysis we will fix these parameters (c s , f2 r o ; z is) to the concordance values of the standard 
model. The only free parameters to be adjusted will then be Q m o and H Q . 

The result of a joint analysis of SNIa, BAO and the position of the first peak of the CMB spectrum is shown 
in Figure [T] In the left panel we show the superposition of the corresponding confidence regions in the (h,Q m o) 
plane, where h = i?o/100(km/s)/Mpc is the dimensionlcss Hubble constant. The central panel shows the same for 
the ACDM model. In the right panel a comparison between the two models is made, after marginalizing over h. 
In this analysis we have used the SDSS SNIa compilation with the MLCS2k2 light-curve fitter [17[- Other current 
compilations make use of Salt or Salt2 fitters which, however, are not model independent. A more complete analysis, 
including the Union2 and Constitution compilations, can be found in jl5| |. In the case of the Constitution sample 
calibrated with MLCS2k2, the results are very similar to those presented here. Note that for the standard model the 
relative matter density is higher than the usually accepted concordance value (£l m o ~ 0.27). This is a characteristic 
of the MLCS2k2 fitter, also present when wc use it to calibrate the Constitution sample. 
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FIG. 1: Left panel: superposition of the confidence regions from SDSS (MLCS2k2) SNIa, BAO and the position of the first 
peak of the CMB in the present model [l5|]. Center panel: the same for the spatially flat ACDM. Right panel: The 
corresponding likelihoods when h is marginalized. 




FIG. 2: Left panel: the 2dFGRS power spectrum in the present model with Q. m o = 0.45 (blue solid line) and for the BBKS 
transfer function with Q. m o = 0.27 (red line). Right panel: superposition of the four tests (SNIa, BAO, CMB and LSS) 
performed with the present model. The blue ellipse corresponds to the concordance region. 



V. THE MATTER POWER SPECTRUM 



The analysis of the matter power spectrum was performed in Ref. [la ], where, for simplicity, baryons were not 
included and the cosmological term was not perturbed. In a subsequent publication a gauge-invariant analysis, 
explicitly considering the presence of late-time non-adiabatic perturbations, has shown that the vacuum perturbations 
are indeed negligible, except for scales near the horizon [l||. Therefore, the result of Ref. [l8| can be considered a 
good approximation. An updated analysis is shown in the left panel of Figure [U together with the 2dFGRS data 



points [20( and, for comparison, the concordance ACDM fitting as provided by the BBKS transfer function [21] . It is 
also shown the 2a confidence region for our fitting. 



VI. JOINT ANALYSIS 



In the right panel of Figure [2] we show the superposition of the la and 2a confidence regions for the four tests 
we are considering here (SNe la, BAO, first peak of CMB and LSS), as well as the concordance ellipses of our joint 
analysis. The concordance value for the matter density is 0.44 < f2 m o < 0.46 (2a). The Hubble parameter is smaller 
than in other estimations, but it has the same value obtained when the same SNIa sample is analyzed for the ACDM 
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TABLE I: Limits to fi m0 (SNe+ CMB + BAO+LSS). 

A(i)CDM ACDM 

Union2 (SALT2) 0.420+° °% 1.063 0.235 ± 0.011 1.027 

SDSS (MLCS2k2) 0.450tg;g^ 0.842 0.260±g;gH 1.231 

Constitution (MLCS2k2[17]) 0.450±g;^| 1.057 0.270 ± 0.013 1.384 

"Error bars stand for 2u 

model (see Fig. [1]). In what concerns the matter density parameter, it is about 10% higher than the flat standard 
model value for the SDSS (MLCS2k2) compilation (0 TO o ~ 0.40). In both models the obtained values of O to q arc 
higher than the usually accepted concordance value for the flat ACDM model, and the same occurs when one uses 
the Constitution sample calibrated with MLCS2k2. A discussion about the present tension between different SNIa 
samples and fitters, in both models, can be found in [l5j . Here we show, in Table I, the best-fit results for f2 m o (with 
h marginalized) with three samples of supernovas: the SDSS and Constitution calibrated with the MLCS2k2 fitter, 
and the Union2 (which is calibrated with SALT2). In Figure [3l on the other hand, the reader can find the evolution 
of the scale factor (left panel) and of the relative energy densities of radiation, matter and dark energy (central and 
right panels) in the two models. 

The concordance value obtained for Q m o is in agreement with estimates based on peculiar velocity measurements 
for galaxy pairs [22[, but only in marginal agreement with dynamical estimations from clusters (2^|. It is important 
to emphasize, however, that such estimations are not very precise and that they involve a different scale. Note also 
that the present result has an approximative character. As discussed in (l8| . the inclusion of baryons in the analysis 
of the power spectrum may lead to a difference of about 10% in the estimation of f2 TO c 

With the concordance values of h and S] m o in hand, we can obtain the age of the Universe, as well as the redshift 
of transition between the decelerated and accelerated phases. They are given, respectively, by 

H t 2 In a2^o (14) 

3(^m0 — 1) 



and 



z T 



'if- 1 



2/3 

-1. (15) 



This leads to H t a = 0.97 and Zy = 0.81, in good agreement with standard predictions and astronomical limits [24j . 
For h « 0.7 we have t « 13.5 Gyr. 

VII. THE CMB ANISOTROPY SPECTRUM 

It is not easy to perform a full analysis of the CMB anisotropy spectrum with alternative models, specially if it 
involves new processes like particle production, for example. Nevertheless, for at least some of the features we do 
not expect significant changes. First of all, radiation and baryons are assumed to be separately conserved (there is 
no photon fluctuations from the vacuum at the tree level, and baryons are too heavy to be produced at late-times). 
Therefore, the CMB temperature, the acoustic horizon at the last scattering and the relative height of the peaks will 
not be affected. 

On the other hand, the redshift of equality between matter and radiation in our model is given by 

*« = ^ (16) 

where we can note again an extra factor f2 m o as compared to the standard ACDM result. By using our concordance 
matter density J7 m o = 0.45, we obtain z eq « 2500, which does not differ too much from the standard value. (This is 
simply because 0.25 = 0.5 2 . Therefore, a relative matter density of around 0.5 in the present model leads to the same 
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FIG. 3: Left panel: The time evolution of the scale factor in the spatially flat ACDM and in the present model. Center 
panel: The density parameters of matter and dark energy in both models. Right panel: The density parameters of matter 
and radiation in both models. In all cases we are using the best-fit values of f2 m o from the SDSS analysis, and the density 
parameters are defined as Qi = pi(z)/3H 2 (z). 




FIG. 4: Left panel: the density contrast as a function of redshift for the Einstein-de Sitter solution (blue), the spatially flat 
ACDM model with Q m = 0.27 (yellow) and the present model with Q m = 0.45 (red), for the same initial condition 6 ~ 10 -5 
for z = 1100. Right panel: the corresponding gravitational potential as a function of z, for k = 0.002 ft/Mpc. 



z eq of the standard model with VL m 0.25). For this same reason, the matter density at the time of equality will be 
approximately the same in both models. [By the way, this result for z eq guarantees that the turnover of the matter 
power spectrum is correctly placed (see Fig. [2] and [18j)]- As the redshift of last scattering is exactly the same 14 1, 
and the recombination does not depend on a particular dark energy model, both last scattering and recombination 
will occur in the matter dominated phase. All these results indicate that the absolute height of the peaks will not 
change either. 2 

The only important consequence of matter production on the CMB spectrum could be owing to the integrated 
Sachs- Wolf effect (ISW). In the right panel of Figure 2] we plot the gravitational potential fc 2 $ = p m a 2 S/2 (5 is the 



2 In reference [SEIl the authors claim, by using an adapted version of CAMB, that a good fitting for the CMB spectrum can be obtained 
with our concordance matter density, Q m o ~ 0.45 (see the left panel of their Figure 3). On the other hand, they claim that this density 
leads to a matter power spectrum in disagreement with LSS data, contrary to what we have obtained in llgll . Our present results, 
however, agree with those obtained previously in reference [lg| | . 
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matter density contrast and k is the comoving wavenumber) as a function of the redshift, for the Einstcin-dc Sitter, 
ACDM and A(i)CDM models, for k = 0.002 h/Mpc. The ISW effect is proportional to the integrated change of the 
gravitational potential. As the right panel of Fig. |4] shows, the slope of the potential of our model differs from the 
corresponding slope of the ACDM model. But also here one may expect this effect to be small. In the ACDM model 
the potential decreases by 24% between the last scattering and the present time, in ou r A( t) model it decreases by 20%. 
Note, however, that Fig. [4] is preliminary since it is based on a Newtonian analysis jl8j . Here, a more quantitative, 
fully rclativistic analysis is necessary. 

On the other hand, as can also be seen in Fig. 2J there is a small difference in the gravitational potential as compared 
to the standard model, a difference that will affect the normalization of the CMB spectrum, since the Sachs- Wolf 
effect is proportional to <£>. But the gravitational potential is proportional to the matter density contrast as well. 
Therefore, if the latter is correctly normalized (by using, for example, the BBKS transfer function [2l|, as was done, 
e.g., in [18[ and in the left panel of Fig. [5]), the CMB spectrum will also have the correct normalization. In other 
words, in order to have the correct normalization today, we have more power in the CMB spectrum at high redshifts 
as compared to the standard case. 

VIII. THE COINCIDENCE PROBLEM 

The coincidence problem within the ACDM model is usually phrased as: why p m and A are of the same order of 
magnitude just at the present epoch? Our model relies on dark-matter particle production at a constant rate. From 
the outset, this rate is a free parameter. By Eq. ((5J it is assumed here to be of the order of the Hubble constant. For 
our model the coincidence problem translates into: why the rate Y is of the order of the present-time Hubble rate Ho 
(and not of the order of the Hubble rate at a different time)? Any other value of Y would not be compatible with the 
observed cosmological dynamics. 

Compared with the ACDM model, the coincidence problem is alleviated, however, in the sense that fi m /f2A oc a~ 3 / 2 
[l9| . whereas in the ACDM model one has fl m /il\ cx a~ 3 . Any decay of the ratio Q. m /Q,^ that is less than a~ 3 can 
be considered as an alleviation of the coincidence problem [26| . 

This difference in the background dynamics also influences the behavior of the perturbations. As the left panel 
of Fig. [4] shows, the present model deviates from an Einstein-de Sitter universe already at higher redshifts than the 
ACDM model. Moreover, while the perturbations within the latter approach a constant value for decreasing z, in our 
model there appears a maximum around the present epoch and perturbations are suppressed for z — > — 1. 

IX. CONCLUSIONS 

The origin and nature of the dark sector is probably the central theoretical problem of modern cosmology. The 
non-zero tiny value of the cosmological term and its coincidence with the matter density are some important faces 
of the question. Relating dark energy with vacuum fluctuations is an old difficulty, but one may expect a potential 
solution if this problem is appropriately considered in curved spacetimc. Several models of dynamic vacuum have been 
suggested along these lines. Models for accelerated expansion based on matter creation in the expanding background 
have previously been proposed, e.g., in |27| . Although a relation between dynamical dark energy and particle 
production is not mandatory, we think it represents an interesting option. 

This paper is based on the assumption that the late-time spacetime expansion can extract non-relativistic particles 
from the quantum vacuum. In the case of a constant creation rate, the backreaction of this process leads to a 
vacuum density linearly proportional to the Hubble parameter, a result previously obtained from estimations of the 
gravitational contribution of the QCD vacuum condensate. The resulting cosmological model is shown to be in 
accordance with the most precise observational tests, namely those which involve SNIa, the anisotropy spectrum of 
the CMB, the LSS distribution and BAO, albeit for higher values of the current matter density than in the standard 
model. Another difference to the ACDM model may be a modified ISW effect. This will be a subject of future 
research. 
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